Stretching an adsorbed polymer globule 
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Using molecular dynamic simulation, we study the stretching of an adsorbed homopolymer in a 
poor solvent with one end held at a distance Ze from the substrate. We measure the vertical force 
/ on the end of the chain as a function of the extension Ze and the substrate interaction energy w. 
The force reaches a plateau value at large extensions. In the strong adsorption limit, we show that 
the plateau value increase linearly in w in good agreement with a theoretical model. In the weak 
adsorption limit, a polymer globule with a layered structure is formed and elastically deformed when 
stretched. In both cases a simple theoretical model permits us to predict the relation between the 
necessary force to fully detach the polymer and its critical extension. 

PACS numbers: 61.41. +e,87.15.Aa,36.20.Ey 



With the new development of single molecule exper- 
iments, the nano-manipulation of individual polymer 
chains and biological macromolecules is becoming a very 
important subject in order to understand their mechan- 
ical properties and characterize the intermolecular inter- 
actions. Forces on the scale of the pico-newton have been 
measured with imposed deformation on the scale of the 
nanometer with atomic force microscope (AFM) 1] or 
optical tweezers pj. Experiments with AFM have been 
conducted with different kinds of biological molecules 
such as DNA Q , Titin 0| and more recently Myosin Q . 
Although experiments on a flexible polymer are more dif- 
ficult to conduct due to its small persistence length, there 
have been several recent studies in both good and poor 
solvent conditions 0, • They show that a stretching 
experiment can give insight on the determination of dif- 
ferent mechanisms of adsorption and on the intermolec- 
ular forces inside the macromolecule. On the theoreti- 
cal side, the stretching of chains in a poor solvent has 
been studied by Halperin et al _8j who predicted that 
at weak extensions the spherical globule deforms into an 
ellipse and that for a finite stretching force there exists 
a sharp first order unwinding transition between a glob- 
ule and a chain. Numerical and theoretical studies have 
confirmed this scenario and have shown that there exists 
a co-existence state between a globule and a stretched 
chain [ilii|ll|ll|ll[ll|ll|ll. As shown in [13 
the critical value of the force saturates toward a finite 
value when the polymerization index N goes to infinity. 
In the presence of an adsorbing substrate force-extension 
relations have been obtained theoretically for an ideal 
gaussian chain [13] , for a poly-electrolyte in good solvent 
[l8| and for a polymer chain in a good solvent [T3I . More 
recently, experimental results have been obtained for a 
stretched polymer chain in a poor solvent adsorbed on 
a substrate [3. In this paper, we consider the problem 
of an adsorbed polymer globule in a poor solvent with 




FIG. 1: Snapshop of the polymer adsorbed on the surface 
with a fixed end (Simulation at lu = 4, T = 0.8, A^ — 100 and 

Ze = 15). 



one end fixed at a distance z^ perpendicular to the ad- 
sorbing substrate. The substrate is characterized by an 
homogeneous interaction energy w with the polymer. We 
mainly investigate the effet of strong adhesion w, how- 
ever we have also perform simulation for intermediate 
adsorption for which the onset of poor solvent effect are 
present. We do not investigate here the weak adsorption 
regime for which most of the chain can evaporate from 
the substrate. Simulation results for the retracting force 
on the fixed end as a function of Ze and w are compared 
with an analytical model. As we stretch the polymer, 
we observe that the force rises until it reaches a plateau 
value which depends on w. 

We use a molecular dynamics simulation technique 
where the equations of motion are integrated with the 




FIG. 2; Graph of /p as a function of w for T = 0.8. Circles are 
simulation results of the plateau value. The line corresponds 
to the theoretical prediction obtained by minimizing F with 
respect to Nc- The inset is a force-extension diagram for 



Verlet algorithm and the temperature fixed with a 
Nose-Hoover thermostat [23 • The interactions between 
monomers consist of two parts : V{r) = Vi(r) -I- V2{r) 
(r is the distance between two monomers). Vi{r) — 
a(r — do)^ + ^(''' ~ ^o)^ is the valence interaction be- 
tween nearest-neighbors (NN) monomers and V2(r) = 
4e((-)^^ — { — )^) the Lennard- Jones interaction between 
non NN monomers which englobes the poor solvent ef- 
fect. Here for simplicity surface tension effects do not 
depend on the polymer configuration, although this de- 
pendence can be taken into account as done in |2ll |. us- 
ing many body forces. Finally, we simulate the pres- 
ence of the atomically flat substrate by adding an at- 
tractive interaction between the wall and the monomers 
: yw{z) = w ((^)^ — (^)^)j where z is the distance to 
the wall. In the following we take e = a = do = 1, 
a = 3000 and b = 10000. The energies and distances 
are respectively expressed in units of e and a and the 
time step is At = 10"'^. For the present study we take 
kbT = 0.8 so that the simulated polymer is under poor 
solvent conditions and still above the freezing tempera- 
ture zV\. We simulate a polymer adsorbed on the sub- 
strate with one end fixed at a distance z^ from the wall 
(Fig. 1). In a first run, starting from a free adsorbed 
globule, a constant force (usually /=10) is applied to the 
end of the chain until it reaches the value Zg. A second 
run, of roughly 10^ molecular dynamics (MD) steps with 
Ze fixed, ensures thermal equilibrium. Finally, during a 
third run of 10^ MD steps, different quantities like the 
time-averaged force on the end of the chain /, the density 
and energy profiles are recorded. As shown on Figure 1, 
the adsorbed globule has a semi-spherical shape due to 
the effect of the poor solvent. 

We consider a polymer composed of A^ = 100 monomers 
for which we perform a set of force-extension measure- 
ments for a broad range of values oi w {2 < w < 80) 
and for Zg varying between 3 and z^, the critical value at 
which the polymer fully detaches from the substrate. For 
all w, we observe that the force rises as we increase z^ 



and saturates to a force plateau value fp for the large Ze 
values. This is illustrated in the inset of Figure 2 where 
a full force-extension curve is plotted. It is then possible 
to plot fp as a function of w (Fig. 2). We observe that 
in the strong adsorption regime, the dependance of the 
plateau in w is linear. For weaker w values, the poor sol- 
vent effects become prominent and the relation between 
fp and w deviates from linearity. 

In Figure 3, we present a force-extension diagram in the 
regime of strong adsorption {w = 60) that presents pro- 
nounced oscillations. For these large values of w, the 
un-stretched part of the polymer is fully squashed on 
the substrate and forms a bi-dimensional (2D) globule. 
These oscillations are due to the discrete character of 
the chain and the peaks correspond to the transfer of 
a monomer from the globule to the chain. As we in- 
crease Ze, f increases rapidly until the next monomer 
can escape from the 2D globule. The stretching tension 
is then released, and / diminishes abruptly. As we still in- 
crease Ze this process repeats as long as there are enough 
monomers adsorbed. Nc being the number of monomers 
in the chain, we furthermore observe a dynamical coex- 
istence between states with Nc and Nc + 1 monomers in 
the chain. This is illustrated in the inset of Fig. 3 where 
we plot the dynamical variation of / together with its 
associated probability density function. For a fixed value 
of Ze — 18, the system can be in two different states with 
A'c = 18 and Nc — 19. As a consequence the distribu- 
tion function of / is bimodal. For smaller values of w 
the morphology of the adsorbed polymer is qualitatively 
different. As shown in Figure 1, in the weak adsorption 
regime the polymer forms a 3D globule partially wetting 
the substrate. As we will discuss below, the top of the 
globule which is in contact with the chain is elastically 
deformed. As a consequence, the necessary force to ex- 
tract a monomer varies continuously and the oscillations 
in the force-extension curve almost disappear (inset of 
Fig. 2). 

We present a phenomenological model in order to un- 
derstand the different results reported above. We assume 
that the fixed force and fixed z ensembles are equiva- 
lent. This is in principles valid, only in the thermody- 
namic limit. We therefore assume that the system is big 
enough to have small fluctuations. As shown in Fig 1, 
the globule-chain system can be described as a stretched 
chain in co-existence with an adsorbed globule. There 
are two main contributions to the free energy F of the 
system : the free energy of the stretched chain Fc and 
the free energy of the adsorbed globule Fg. Let Ng and 
Nc be respectively the number of monomers in the ad- 
sorbed globule and in the chain, so that N = Nc + Ng. 
The major contribution to Fc comes from the reduction 
of the fluctuations of the stretched chain |23| 

Fc = -NchT log U^j^ sinh if do /kbT)] + fzc , (1) 

where the potential energy of the chain is neglected. The 
dominant part of Fg comes from the interaction energy 




FIG. 3: Force-extension diagram in the strong adsorption 
limit for T = 0.8 and w = 60. Circles are simulation re- 
sults and the line is the theoretical prediction presented be- 
low. Inset: Ze — 18, time variation of the force together with 
its associated probability density function illustrating the dy- 
namical co-existence between the two states with Nc — 18 
and Nc = 19 monomers. 



of the monomers with the substrate and the interaction 
between the monomers. We can neglect here the en- 
tropic contribution to Fg |23|. In the regime of strong 
adsorption, the interaction between monomers within 
the 2D squashed globule can be neglected relative to 
the substrate interaction. For large w values we have 



Fn 



-NgW„ 



where Wmin is the substrate potential 
minimum ( Wmin — O.SSw). We write the total free en- 
ergy F = Fc + Fg where z^ is fixed and / is an implicit 
function of Ze and N^. Here Zg is given by the non-linear 
Langevin relation |J22J : 



Ze = Ncdo[coth{ 






kbT 



(2) 



We express the free energy F as a function of the un- 
known Nc by inverting numerically equation ^ for /. 
We then minimize numerically the free energy F with re- 
spect to Nc by keeping w and Ze fixed. For a large value 
of Ze, we substitute in equation (O the obtained value of 
Nc to get the plateau force fp. The simulated fp values 
presented in Figure 2 are in good agreement with our 
theoretical model for large w. In the high w limit and 
for r = 0, we expect fp to scale linearly with w so that 
Ze = Ncda. In this case the entropic contribution can be 
neglected relative to the free energy of the globule energy 
so that fp oc 0.38i(j in relative good agreement with the 
simulations {fp oc OAOw for T = 0.8). In the weak ad- 
sorption limit, the model underestimates fp since the in- 
teraction between the monomers has not been included in 
Fg . In this case one should add to Fg monomer-monomer 
interactions of the wetting globule together with surface 
tension corrections. The surface tension correction gives 
a contribution of the order of j{N — Nc)'^^^ to Fg, where 
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FIG. 4: Force at the detachment point z^ for different values 
of w and N. Diamonds: A/" = 100 w = 3, 4, 6, 8, 10, 20, 40, 60 
from left to right. Circles: A^ = 200 w = 4,8,10,20,40. 
Square: A'^ = 400 w = A. Full lines are theoretical predictions. 



7 is a phenomenological surface tension parameters. This 
contribution plays only a minor role. We find by perform- 
ing the same minimization as above that there is a very 
slight decrease of the plateau value as Ze increased, this 
is also confirmed by our MD simulation. A good agree- 
ment between simulation and theory is also found for the 
extension-curve of Fig. 3 obtained for various values of Ze 
using the same procedure as above. The smoothing out 
of the peaks is due to the dynamical coexistence between 
states with Nc and Nc + 1 monomers in the chain. This 
dynamical coexistence is induced by the thermal fluctu- 
ations and gives rise to important time variations of /. 
As we increase the value of Ze , we assume that the poly- 
mer fully detaches from the substrate when Nc = N . 
Using the relation between /, Zc and Nc given in equa- 
tion (|2Jl, we obtain the relation between the force / and 
the critical extension z^ at which the polymer detaches 
from the substrate. We plot in Fig. 4 the theoretical 
curve f{ze) and the simulation results for different val- 
ues of w and different system sizes. Here again a good 
agreement is found between theory and simulation for 
the large w values. The agreement is slightly less satis- 
factory for weaker w values. In this case, we observe that 
the detachment occurs suddenly when a finite number of 
monomers are still adsorbed on the substrate (TVc < A^). 
This is analogous to the first-order unwinding transition 
observed on a globule under tension and discussed in the 
introduction. 

We finally investigate the globule structure in the regime 
of weak adsorption in order to quantify its elastic de- 
formation upon stretching. We represent in Fig. 5 the 
monomer density p{z) and energy E{z) profiles along z 
for different values of Ze- The density p{z) presents de- 
caying density oscillations similar to the one observed for 
a liquid in the presence of a solid substrate [2J| . The pres- 
ence of the density oscillations are due to the poor solvent 
effect. We have checked that in the good solvent case there 
are no density oscillations. As shown in Fig. 5 (top) the 
overall shape of the density profile in the globule is in- 




FIG. 5: top:Ai' — 100, density profiles p{z) showing layering 
for a free globule and globules with Zg = 15 and 20 (w = 8). 
Bottom: energy profiles E{z) for the same parameters. 



variant, up to a normalization factor, as we stretch the 
chain. The layered structure can therefore co-exist with 
the stretched chain. As the monomers are extracted, the 
globule shape remains qualitatively the same, indepen- 
dent of the size and of the tension as long as the number 
of monomer in the globule is large enough. In Fig. 5 
(bottom), we plot the energy density profiles E{z) where 
E(z) is the sum of the valence and of the Lennard- Jones 
interaction energies. We can see from E{z) plotted for 
different values of Ze that there exist three different re- 
gions. In the first region (z < 3) mostly corresponding to 
the first adsorbed layers, the energy is independent of Zg. 



This is consistent with the above discussion on the negli- 
gible globule shape variation upon stretching. A second 
region, roughly comprised between z — 3 and z = 8, is 
formed between the globule and the chain. In this match- 
ing region the top of the globule is elastically deformed 
and pushed up by the tension. The overall vertical me- 
chanical equilibrium is ensured by the elevation of the 
top of the globule. The small difference in the energy 
profiles computed for Zg = 15 and Ze = 20 in the second 
region also demonstrates that the elastic deformation is 
independent of the stretched chain length. Finally the 
third region (z > 8), corresponds to the stretched chain 
with a very small potential energy, confirming our main 
theoretical hypothesis. 

The picture we have presented above holds also for 
larger value of N. Simulations which are not presented 
here have shown that fp depends only weakly in N and 
finally saturates. To summarize, our simulations have 
shown that there is a plateau value for the force that 
depends on the substrate interaction w. In the strong 
adsorption limit the un-stretched part of the polymer 
forms a 2D globule. A dynamical coexistence between 
states with Nc and Nc + 1 monomers in the chain has 
been observed. For weaker w values the adsorbed poly- 
mer forms a 3D globule wetting the substrate and elas- 
tically deformed upon stretching. A simple theoretical 
model of the force extension relation is presented and 
leads to a good agreement with simulation results. The 
infiucncc of the chain persistance length, of the tempera- 
ture and of finite size effects will be presented elsewhere. 
We hope that this work will motivate new experiments on 
the stretching of long molecules and bio-molecules. The 
weak adsorption regime will be investigated elsewhere. 
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